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ABSTRACT 


Theoretical analyses of flow with rotating passages 
or turbo-machinery have become necessary for the proper 
design of turbo-pump elements in liquid fuel boosters and 
power conversion units, One such theory is presented, and 
a numerical solution derived, This theory develops a meth- 
od for the analysis of steady, inviscid, adiabatic flow 
through arbitrary rotors, A detailed analysis in a merid- 
ional plane is given, assuming axial symmetry, A simplified 
approach to the blade to blade solution is also presented, 
The merits of these theories are compared with other pro= 
posed methods, The inverse, or design, approach is con- 
sidered, and found to be unnecessary, 

A numerical solution for incompressible flow is derived 
and applied to the flow solution in the impeller of a mixed 
flow compressor with backwards-bent blades cf arbitrary 
shape, Meridional streamlines and relative velocity dis- 
tributions are progressively calculated on a CDC 1604 com- 
puter, using FORTRAN program language, Data are measured 
from a detailed presentation of the blade shape in a merid-= 
ional plane, Blade to blade relative velocity distributions 
are calculated from the meridional plane analysis, 

It is concluded that the results completely define the 
flow and are sufficiently accurate for engineering applica- 
tions, Validation is based upon the reliability of the the- 
ory, and upon comparisons with results of other methods, 
Extensions of the scope of this approach are recommended, 
which include the compressible solution and the solution of 
flows in unbladed passages, 

The writer wished to express his appreciation for the 
assistance and encouragement given him by Professor M, H, 
Vavra of the U. S. Naval Postgraduate School, His guidance 
in the development of this solution was invaluable, 
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A METHOD FOR THREE - DIMENSIONAL 
FLOW ANALYSIS IN A ROTOR 


USING A HIGH SPEED DIGITAL COMPUTER 
INTRODUCTION 


There is an increasing demand for compact, high speed, 
high power output turbomachines for use in liquid fuel roc- 
ket engines and power conversion units, These units are 
being designed to work with such media as cryogenics and 
liquid metals, The proper design of such machinery largely 
depends upon the accuracy of the theoretical analysis, used 
prior to fabrication, The method used should reasonably 
predict aerodynamic forces and tendencies toward flow sep- 
aration and cavitation, 

The more simplified approaches, using conditions ahead 
of and after the rotor, ignore conditions within the rotor 
passage itself, and therefore, are unsatisfactory, Poten- 
tial flow analyses are useful for stationary cascades, but 
cannot be accurately applied to flows in rotating passages, 
The need for more precise methods, required for the proper 
analysis of modern rotor designs, was anticipated, Several 
three-dimensional theories have been postulated, In gen-= 
eral, these theories are quite complex and are difficult 
to apply. 

One such theory was developed by the NACA's Lewis 
Laboratory in the early 1950's, The three dimensional 


problem was reduced to an iterative process between two- 
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dimensional solutions on hub to tip, and blade to blade 
stream surfaces, Ref. 1 presents an analysis on a merid- 
ional, or hub to tip, surface, A similar approach is de- 
veloped in Ref. 2. A blade to blade analysis is described 
in Ref, 3. These treatments were combined into a general 
theory in Ref, 4, Subsequent efforts were directed towards 
applications of this theory to specific examples, The com- 
plex nature of the theory reguired the introduction of cer- 
tain simplifying assumptions, The development of the high 
speed digital computer made solutions more practical, as 
demonstrated in Refs, 5 through 7, 

Prior to the NACA's work, a method of solution ina 
meridional plane, similar to that of Ref. 1, was developed 
by Meyer in Ref. 8. This method uses an iterative, graph- 
ical process to solve two simultaneous, linear differential 
equations along the characteristics of the equations, The 
equations are derived from the Eulerian equations of motion, 
and the characteristics are defined by the geometry of the 
flow channel, In Ref, 9, Vavra reorganized Meyer's scheme 
into a more general theory, reducing the equations of mo-]= 
tion to one non-linear equation which is also solved along 
characteristics, In addition, a simplified blade-to-blade 
analysis was developed, 

The purpose of this thesis is to transform the theo- 
retical development of Ref, 9 into a method of solution 
for arbitrary rotors, using the CDC 1604 digital computer, 


A completely contained computer solution was presumed to 
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be quite long, It was anticipated that the amount of in- 
put data, required for acceptable accuracy, would exceed 
computer capacity and that long preparation and computer 
run time would be required, It was decided to use a ser= 
ies of short computations, which could be repeated as often, 
or at such intervals, as required to obtain the desired 
accuracy. 

The theoretical development may be applied to both 
the design and inverse problems, The major portion of this 
thesis is devoted to the latter application; the analysis 
of flow in a given rotor, A design attempt is made to pro= 
vide a physical model for the evaluation of the analytical 
methods, This attempt is not completed because of time 
considerations, and an actual rotor is used for the model, 
This impeller is a part of a compressor test rig located 
at the U, S, Naval Postgraduate School, Details of this 
machine are enumerated in Ref, 10, The rotor is of the 
mixed flow type, with non=-radial blades, The blades have 
a constant thickness, and are of the deloaded type, with 
a reversal in curvature over the rear section, This con- 
figuration presents an arbitrary design which does not 
conform to standard impeller types, 

The methods developed in this thesis are applied to 
this model, using simplified initial conditions, which do 
not reflect actual operating conditions, but are within 
the compressor's operating range, No correlation is made 


between theoretical and test results, however theoretical 
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results are presented for validation, Both the theory 
and the results of this analysis are compared with those 
of Refs, 1 through 7, to substantiate the applicability 


of this method, 





II, THEORY 


A complete derivation of the theory is presented in 
Ref, 9. This section contains those equations and devel- 
opments that are considered necessary for clarity and 
continuity, in following the subsequent transformation 


from theory to numerical methods, 


A, Meridional Plane Analysis 

1. The following assumptions are made to establish a 
model for the hub to tip analysis, A more thorough treat- 
ment of these assumptions is conducted in Section IV, 

a. The rotor cascade contains an infinite number of 
infinitely thin blades, Thus, stream surfaces and 
fluid motion are axisymmetric, 

b, Flow is inviscid, steady, and isentropic, 

¢c, Entropy changes due to discontinuities at the lead- 
ing and trailing edges are acknowledged but ignored, 

a, Flow is incompressible, 

2. The solution of the equation of motion is based upon 
complete definition of the geometry of the blade surface, 
A set of meridional streamlines, m, are assumed and the 
normals, n, are drawn, establishing an orthogonal, axisym- 
metric coordinate system, This system is corrected by 
successive approximations. The blade surface is repre- 
sented by the circular projection on the meridional plane 
of the lines of intersection of the blade surface with 


planes O=constant, Ois the angle measured in the 
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peripheral direction, The © constant planes are planes 
extending in a radial direction from the axis of rotation, 
perpendicular to that axis, The peripheral angle, © , is 
referenced to an arbitrary point on the blade surface, 
usually on the leading edge, The © =constant curves are 
shown in Fig. 1. Assumed coordinate systems are illus- 
trated in Fig. 4. 

The angles defined by this family of curves, at an 
arbitrary point, P, are shown in Figs, 2 and 3, The angle 
& is the inclination of a © =constant curve from the nor- 
mal, The angle A is the inclination of a meridional 
streamline from the axial direction, Z ., The angle (3 
is the flow angle of the relative velocity, W, on the 


stream surface, 
as a9 
tan B=R at (1) 


The deviation of a blade section, € , from the radial 


tan € = tan @ sin (A- 5) (2) 
Cos 


Thus, the system of streamlines and © =constant curves 


direction, R, is: 


are sufficient to completely define the blade shape, 
The relative velocity on a streamline may be ex- 


pressed vectorially as: 


W=O Wy ton@ + Yo Wy (3) 
wheres: WU = unit vector in peripheral direction 
mM = unit vector in stream direction 


The flow problem is reduced to a solution for the merid- 


ional component of the relative velocity, W.. 
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3. The Eulerian equation of motion is: 


a+ 0H = Vx(VxV) + Tys+5 (4) 


For steady, isentropic, relative flows: 


Vp = Wx (yx«W +2 O)+5 (5) 
where, 
a Ww? 132 
He = ht Ww" — eR (6) 


For inviscid flow, the friction force, f, is zero, How- 
ever, the effect of infinitesimal pressure changes across 
an infinite number of blades is accounted for by intro- 


ducing the blade force, F,, which is normal to a blade 


element, Thus: 
VHp = Wx (Wx W426) +4 Fa (7) 


Eq. (7) is reduced to its scalar compenents, using: 


of ae (8a) 


a N 
Fa = Fy (a - WH tan@ + NW tan ton 5 ) (8b) 
wheres: A = unit vector in normal direction 


The G component yields: 


RFu = Wy 2 (R Woy + OR") oo R*) 


ae (9a) 


The it component, with Eq, (9a) gives: 


OHe = 
am 0 (9b) 
This relation only holds within the rotor for the assumed 
theoretical model, Flows outside the rotor are unaffected 


by such a model, 
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For all but design conditions, there is a flow dis- 
continuity at the leading edge, since flows ahead of the 
rotor will not meet the rotor at blade entry angles, It 
has been assumed that entropy changes caused by these dis- 
continuities are neglected by ignoring the changes in total 
enthalpy across the leading edge, Conditions ahead of the 


rotor are derived from Eq. (4) for steady, isentropic flow, 
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A relation between conditions on a streamline ahead of the 
rotor and within the rotor is established as: 


He = ota OOK, (\Nuc + coR,)\ (11) 


where e and L denote entrance and leading edge stations, 


The change of Hp along a normal is: 


he = OKe OMe —~ ow O(KRL Wy + W Re ) db 
an Me ON Die aM (12) 
= K(n) 


The N component of Eq. (7) becomes: 


Wy OWNrn + Win Reg + Wu OCRWay + &R*) 
on BR an 
(9c) 


+ Wm O(CRWy+ UR") tan @tan oO = Kin) 
KF om 


where Kn is the curvature of the meridional streamline, 


From Fig, 3: 
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Eq. (9c) is reduced to a relation in wa! 
ve awe 2 
oe (1++an@) + em (tan‘*a tans) 
Z O(Rta o(R tan GB) 
+Wrn (2Rm+ 2 tone | oiktan®) + tan § o(R tan 8) ) (14) 


+ Wr (4 os tan@ Q& +tand QR | \ 
= 2K(n) 


The partial derivatives in brackets are reduced to a total 


derivative along a © =constant curve by: 


YL Se oC) Ol ) 
che Sin § Ss + coss oe (15) 


Eq. (14) becomes: 


2 2 if 2 2 
= + sin? tend 2Win + Wm (2Ry Cos G 


a sin (28) d (R tan Q) We cdf 2 Sian aR \ (16) 
Q )* ™ (eee 


mCGsam ad al 
=P Cos (a Kn) 


Eq. (16) is reduced to an ordinary differential equation 
along a characteristic, whose tangents at any point satisfy 


the relations: 
ae 
tan ¥ = sin @tan 6 (17) 


According to Ref. 8, these characteristics are regarded as 
circular projections on a meridional plane of unique spa- 
cial lines lying on the blade surface, Introducing Eq. 
(17) into Bq, (16): 
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where x denotes distance along a characteristic, ands: 


ie = cos {2 Ryn Cos + Sin (2A) gee (19a) 


Rcosd 
= Sin (2@) dR (19b) 
Yo = 2cos ¥ wt 4 
Ya = EOS: oh, or ane KC) (136) 


4, The value of K(n) is dependent upon conditions ahead 
of the rotor, as shown by Eq, (12). The normal component 
of the equation of motion for steady, isentropic flows 
ahead of the rotor blading is: 


Vin SEP + Vin Ren + Ye ORY BE = (20) 


Eq. (10) shows that for isentropic conditions at the en- 
trance, the total enthalpy and the product RV, are constant 
along a meridional streamline, This condition is proved 


from the 0 and ™m components of Eq, (4). Thus: 


2 
Nin in X; + Xe = (6 (21) 


x, =2Ry, (22a) 
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The solution of Eq. (21) is: 


~}X,dn X,dyn (23) 
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This solution is solved for given entrance conditions 
and the volumetric flow rate, Q, by iterating the hub 


velocity, V_.,, until the flow rate is satisfied by; 


mH 
r 
Q=2r { RVm tana Aan (24) 
At the leading edge: Bs 


Thus, K(n) is solved for given entrance conditions, 

5. The remaining terms of Eqs. (19) are determined from 
blade geometry and the orientation of the meridional 
streamlines within the rotor, These streamlines must be 
assumed and corrected by successive approximations, The 
location of the streamlines at the leading edge are de- 
termined from entrance conditions, Eq. (18) is solved 


along a characteristic until the flow rate is satisfied by: 


Kop 
Qs an | R Wm Cos 6 dX oll 


o 

B, Blade to Blade Analysis 
1, It is assumed that the relative velocity varies line- 
arly across the blade channel, This assumption is taken 
from thin airfoil theory, where blade profiles are replaced 
by bound vortices, The relative velocity from the axisym- 
metric solution is considered to be the mean velocity 
along the periphery at any point in the meridional plane, 
These assumptions lead to: 

Wp = W- AW/ (27a) 

We = \V + AW (27b) 
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where p and s refer to the pressure and suction sides 
of the blade, 
2. The difference between relative velocities across the 


blade is related to the static pressure difference by: 


P+ Cp \Ww* = Ccons+ant (28) 
Thus s 
- &P 
AW = 20 \N (29) 


The pressure difference across a blade element, which 
is projected onto a meridional plane, is related to the 
moment about the axis of rotation, exerted by the flow on 
that element, by: 

WA eo 
AM = 2(RApdndm) (30) 


This moment is also derived from the momentum theorem and 


Eq. (9a): 
AM = - BIR Aedndme (RW, + wR?) (31) 
- om 


The minus sign indicates that the moment oppose the rota- 
tion of the rotor, The arc, AO , is expressed in terms 
of an equivalent blade thickness, At’ , which is measured 


perpendicular to the meridional plane, 


6 32 
At =~ 4t  ~/j+ sin’ Cos S —_ 
Cos (3 


(33) 
Te Pd Zam At 
N R 
where: N = number of blades 
At = blade thickness 


BZ 





Equating Eqs. (30) and (31) and combining with Eqs, (29) 
and (33): 
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III. NUMERICAL SOLUTIONS 


A, Design or Inverse Solution 
1. The application of the preceding theory to impeller 
design was investigated in an attempt to produce a physi- 
cal model for the direct solution, or flow analysis, Over- 
all dimensions of an axial entry, mixed flow pump impeller 
were provided by Professor Vavra, which included: 

a, Hub and tip profiles of meridional contour 

b, Leading and trailing edge contours in a merid- 

ional plane 

¢c. Blade entry and exit angles 
The problem was reduced to that of determining blade shapes, 
which would not only satisfy these boundary conditions, but 
would conform to practical structural and aerodynamic limi- 
tations, 
2. A drawing was made of the meridional plane anda series 
of approximate normals were constructed, A meridional 
streamline was constructed, dividing the rotor annulus in 
natt. 

ae a (R.+ Em malic (35a) 
Iterations were made on Ra along a normal until; 
Ane (Ry + Ry) M1, (35b) 

The slope of this streamline, \ , was measured, plotted, 
smoothed, and the streamline reconstructed, Normals were 
corrected, 

The distribution of the flow angle, (3 , along the 
mean streamline was assumed in the form: 
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= A+ Brn* (36) 


A and B were evaluated from initial conditions at the lead- 
ing and trailing edges. Values of © were calculated by 


numerical integrations (Ref. 11), From Eq. (1): 


i 


©. 


L 


2 
ein + f tan (3 alm (37) 
FR 


L~| 
where: ) 


L =O 
The slope of the ©=constant curves at any point on 
a streamline may be calculated by Eq. (2), if € is known, 
The choice of an € distribution is governed bys 
a, Maximum € is limited by structural considerations 
b, € distribution must be compatible with fabrica- 
tion tooling practices, An example is given in 
Ref. 9, 
The design problem could have been simplified by using 
radial blades, where € = 0, However, this is contrary 
to the intent of this thesis in presenting methods applic- 
able to arbitrary designs. Therefore, several © distri- 
butions were assumed, and the angles, 8S , computed on 
the mean streamline, 
3, The problem of analytically generating a complete 
family of © =constant curves was as yet unsolved, The 
(3 and € distributions on streamlines other than the mean 


are not independent of the mean distributions, Arbitrary 


choice of these distributions would, in all probability, 
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generate impractical blade shapes, These distributions 
should be assumed as: 

B= aCm,n) 

é€=€E(M yn) 
It appeared that the design problem had become more complex 
than originally anticipated, It was decided that a complete 
treatment would detract from the main objective of this 


thesis, Therefore, the design analysis was discontinued. 


B, Meridional Plane Analysis 
1, The approach used to solve Eq. (18) along a chosen nunm- 
ber of characteristics is, briefly: 

a. Assume meridional streamlines in the vicinity of 
the leading edge which divide the flow into ap- 
proximately equal increments, 

b, Construct normals, thus establishing a coordinate 
grid system, 

c, Measure the necessary blade physical characteris- 
tics at each grid point, 

d, Generate a characteristic curve within the grid 
system, 

e, Compute the coefficients, Y, of Eqs. (19). 

£, Solve Eg. (18) for vee at the intersections of the 
characteristic with the streamlines, iterating 
Wall until the flow rate is satisfied, 

g. Correct these intersections until the streamlines 


divide the flow rate into the prescribed increments, 


h, Recompute ee at the new intersections, 
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i, Project the corrected streamlines farther into 
the rotor channel, and repeat the process for a 
new characteristic, 
Thus, the meridional streamlines are generated from lead- 
ing to trailing edge, and tne distribution of the relative 
velocity on a meridional plane of the rotor is computed, 
Steps a, b, ¢c, and i are solved by the computer, The re- 
maining steps are carried out by graphical means, 
2, There are two categories of initial data required, 
The first is a complete physical description of the rotor, 
Finally, flow conditions ahead of the rotor must be pre- 
scribed, Normally, sufficient data are available from 
drawings of the impeller to construct a meridional plane, 
The © =constant curves are best produced by orthogonal 
projections from a three-view drawing, The meridional 
plane for the impeller used in this analysis is shown in 
Fig. 1. 

The flow conditions at an entrance station must either 
be prescribed from an analysis of the machine installation, 
or assumed, The specific data required are: 

a. Thermodynamic data for the fluid 
b, Velocity distribution 
c¢, Flow rate 
ad, Inlet channel contours 
e. Impeller RPM 


Conditions for this analysis are simplified by assuming: 


Ly 





N= NC = comet 


mL 


DHe _ 
Pays 


This flow is somewhat impractical since it does not 


O 


reflect actual velocity and energy distributions imposed 
by intake ducting or guide vanes, However, these distri- 
butions do present a definite off-design condition at the 
leading edge, Design RPM, and a flow rate compatible with 
actual machine operating conditions, are used, Initial 
data are listed in Table I, The meridional velocity at 
the leading edge was computed by: 

View = = (38) 
A, was computed along the leading edge by Eq. (35b). 
3, Streamlines are extended from initial points on the 
leading edge, such that the flow rate is divided into egual 
increments, Eight divisions were used for this analysis, 
as shown in Fig. 4. 

A grid system is constructed in the vicinity of the 
leading edge. A point on the hub streamline is chosen as 
the starting point of the first characteristic, This point 
is determined By, predicting the approximate alignment of 
the characteristic, Eq. (17) indicates that the slope of 
the characteristic, XY , is proportional to, but less than, 
the slope of the © =constant curve which originates from 
the starting point, % also has the same sign as § ., 
These guidelines fix the characteristic within a desired 
region, These considerations also assist in estimating 
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the proper grid width, A "starting normal" is constructed 
from the starting point on the hub to the tip streamline, 
A second normal is constructed adjacent to the first, in 
the predicted direction of the characteristic, The dis- 
tance between the two normals on each streamline is used 
to fix the remaining grid points, These points do not 
coincide with true normals, however equal grid spacing on 
streamlines is used to good advantage in subsequent calcu- 
lations, 

The values of A , R, and Z are measured at each grid 
point, The values of © were read at the intersections of 
the © =constant curves and the streamlines, plotted, and 
© at the grid points found by interpolation, The angles 

O are determined in a similar manner, The deviations of 
the © =constant curves from the radial direction, 6-nr , 
are measured, \ is added to the interpolated results to 
obtain §. 

4. Each point on a characteristic must satisfy Eq. (17), 
therefore, the angles (3 and 6 must be calculated for any 
point on a streamline, (> is defined by Eq. (1). The 


th streamline is 


derivative of © at a grid point on the j 
approximated by five-point difference formulas from Ref. 12, 


Forward Differences 


dO, = -250; + 48 O14) - 36 Oi424+169;, 53 Ory 398) 
ia ‘2 AN; 


ro 





Central Differences 


(39b) 





Backward Differences 


ASi - Z2Qj-4 16 Oj-3 + 3CSi-72 — 48 On s*235Oy 39.) 


JS M5 I2ZON; 


n streamline, The de- 


AN, is the grid spacing on the ae 

rivative of © at points other than grid points, are cal- 

culated by interpolating linearly between adjacent normals, 
oq =-S2i 4 (Sin = =. (40) 

« ™m 4 AM 5 dm, }) ANy 

GN is the distance from the i°" grid point to the point P. 

A consequence of Eqs, (39) and (40) is that the minimum 

width of the coordinate lattice is six grid points, 

The values of R, \ , and © ata point p™ are also 
calculated by linear interpolation, Q and % are calcu- 
lated by Eqs. (1) and (17). 

5. A characteristic is approximated by a polygon with each 
side terminated by adjacent streamlines, Each intersection, 
ae is located from information calculated at the preceed- 
ing intersection, In this derivation, the lines between 
grid points are assumed to be straight, The derivation is 
shown graphically in Fig, 5, 

Assume that P, is established, and the angle ee 
known, A straight line from ee to P' on the next stream- 


line at the angle “Y . The distance, @N' is: 
AN = aAMtanY% (41) 
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wheres 


aR 
aM = 28 (42a) 
oz $ al 
ee 
OM = 1s (42b) 





Eq. (42b) is used for \ greater than 45 degrees, The 
4 

angle W is calculated by the methods of paragraph 4, 

A line is extended from Po to P'® at this angle, The 


* 
point P is established by: 


P* is the reference point for the extension of tne char- 
acteristic to the next streamline, The polygonal approx- 
imation of the completed characteristic is smoothed through 
the points P*, (Fig, 4) 

There are a number of considerations that must be 
accounted for in generalizing this method for computer pro- 
gramming, dN* originates at the starting normal, The 
Sign of A@N' is the same as that of Y . The sign of 5N 
depends on the difference between XY and ia The sign 
of dN* may be plus or minus, and its magnitude may be great- 
er or less than that of dN’, The value of dN in Eq, (40) 
for a positive dN* is: m 

dN=dN (44a) 
while, for negative dN*; 
aN= AN-dN* (440) 


6. The coefficients, Y, of Eq. (18) contain four elements 


which are, as yet, undetermined, 
a. Ryn 


Za 








b, g (Rtan @) 


an 
c. dR/dQ 


aeri< (n\ 
The curvature, k, of a given curve is extremely difficult 


to calculate, If the curve can be expressed analytically: 


fi 


a 
[\ + (R/)?\7/2 


where R' and R" are the first and second derivatives, 


(45) 


respectively, of the function R = £(Z), 

An attempt was made to determine the equations of 
the hub and tip contours in cartesian coordinates: R = £(2Z). 
Two polynomial approximations of these curves were made 
using CDC Cooperative Library routines, The coordinates 
of the curves at each quarter inch along the axis were in- 
troduced as data, Only high order polynomials fit these 
data points properly, The derivatives of these equations 
were difficult to obtain witnout frequent interruptions in 
computer operation, When calculated, the second deriva- 
tives reflected the sinuous nature of polynomial approxi- 
mations, resulting in inaccurate curvatures which varied 
in sign along a given curve, These methods were considered 
to be unacceptable for this study. 

The derivatives of the hub and tip contours were cal- 
culated with various finite difference equations, (Ref, 12) 
The three and four point calculations were inaccurate, The 
higher order solutions resulted in erratic second deriva- 
tives, It became apparent that all methods were very sen- 
sitive to the accuracy of the data, It was found that 


Zz 





inaccuracies in the third decimal place of radius data 
were sufficient to introduce gross errors in the second 
derivative, 

A semi-graphical approach was used to determine the 
derivatives, The first derivative was calculated by 
finite differences, and the results plotted and smoothed, 
This process was repeated to calculate the second deriva-~ 
tive, using the plot of the first derivative for input, 
Data taken from the two plots were used to calculate the 
curvature in Eq. (45). Hub and tip curvatures are plot- 
ted in Fig, 6, 

This procedure is obviously incompatible with con- 
tinuous computer operations, The computation of the 
streamline curvatures, Kn? at points P*, would be quite 
tedious using this method, Therefore, Kn is approximated 
within the flow passage by linearly interpolating between 
hub and tip, using the data of Fig, 6, Kn is assumed to 
vary in equal increments between streamlines: 


R Mj (46) 


Mr 
Linear interpolation, similar to Eq. (40), is used to cal- 


in Thad Me y= eee 


culate Kn between grid points on a streamline, 
The following derivation defines =¢- (R tan ) in 


terms of known quantities, From Eq. (15): 


d(Rtanga) . «; o(Rtan®) OlR tan) 
ge eB sin § a + Cos & —— (> a 


= ‘Ot ‘ot oF AR 
=k (sin8 Sen? + cos § Zon) + Tan (3 Sy 
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f= 


Prom Fig, 2: ar 
ay ~ COs (G6-A) =D, (48) 
From Eq. (1): 


oieten() =R 22.5 08 OK 





am am 
o(Rt+an @) = Pp O'S ae oR 
on ‘dmon OM ON 
d(Rtan@B) . ~2 cl?6 c\2O 4 
na oo = R (sin eae = COS el Zi LONE 


Grouping terms of Eq. (19a): 





d(Rtan@) — 178, do 
Peas AQ = 5 (tana $8 2 Sendn 
(49 ) 


R cos5 
The term a*@ /am? at a point P* is computed by as- 
suming a parabolic distribution of d@/dm between grid 


points, 





o = a+ bm -+ cm (50) 
+ a8 3 x _ cen G 
_ aA*s *< 
Fig, 7 will assist in following the derivation for b 
™; 58.) = (a4 brn, + Cm?) i$ (52a) 
As a 2 & (52b) 
Pow > (a4 Sra + cm; ) ni} 
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Subtracting iq, Bon from atts 


= me L(am (a&)~|» vn? | (a@ \a (53) 


mM, mF - Lome 


m, and m5 are defined in Fig, 8, 


2 
The function, S , is calculated by differentiat- 
ao 


ing aa along the normal nearest the point P*, A para- 
bolic distribution of — over three streamlines along the 
normal is assumed, The coefficient, b, is derived in the 
same manner as that described in the preceding paragraph, 
In this derivation the two succeeding grid points from the 
origin are used, instead of adjacent points, Calculations 
for points on the last two streamlines are made by using 
the two preceding grid points, 


Forward Differentiation 


5 -_ne WGR),= -S er ni? L(Grn)e ae, | (54a) 


agers: oat ean 
Backward Differentiation 


n2\(Sa)., - ( $2).|- n2\(3 eis -2 (dan) j|(54b) 


ie rales re N12 Male 
The increments n, are functions of the grid spacing, dM, 





ye 


The values of K(n) are required for the calculation 


of Y5. Initial conditions have simplified Eq, (12). With 
Eq. (29) 
w = 1-RPM (56) 
30 
= ~ aR a sm 
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where © is the slope of the leading edge, (8, is cal- 

culated, using the forward difference formula of Eq. (54a) 
in conjunction with Eq. (1). The terms in parentheses in 
Eq. (55) are grouped, and the derivative calculated using 


the parabolic distribution methods of Eqs, (54). 


F(LY = RLV EON BL + Re (58) 
In summary, the coefficients, Y, are reduced to: 
Vine Coe aa 2a cox + SIN (2 nes D2| (59a) 
Yo= 2009 % sim (213) ee (59b) 
CEs 
Yan= 2cosy cos @ Kin) (59c) 


7. The equation of motion, Eq, (18), is solved by succes- 
Sive approximations, until the distribution of Wa along 
the characteristic satisfies the flow rate of Eq. (26). 
The steps in this solution are: 


Pe Assume W 
e Mm 





H 2 
b. Solve Eq. (18) for 3a mn 
cl x 
c, Let: 2 AWN inv 
Wm2 = Wma dx oe (60) 
- aM, 61 
aN] COS My (oy 
GQ, Solve Eq. (18) for 2 
cd Wine 
€., Lect: J x Ne 
ee. dW? W wy 2 
Wren = Wine +f ( yh + SYS \ axe - 


£f. Repeat for Wad! 


The variables under the integral sign in iq. (26) are 


using W_. in (a). 


grouped, 
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E(x) = RW, Cos 


Q = an E(x) dx (63) 


4 
Numerical integrations are performed over two adjacent 


streamlines at a time, and the results are successively 

summed, A parabolic distribution is assumed for E(x), 
E(x\ = a+ bx+C x* (64) 

The coefficients are derived in the same manner as that 


used for the development for the difference equations, 


a= = (65a) 

b — %2 (E,-€.)- x? (Ex-E) (65b) 
Khe — Re XP 

c= &K(E-E)-%, (E2- E.) (65c) 
Ki2X%2- Ke X, 


wheres: 
Ao JX, 


X2 =d\,+ d Xz 
The integral of E(x) is: 


R+2 
Z 4 “2 yt 
E(x = A(X, ,.7 Xk) * Bee me) a 
R 3 3 
c 
+ $ (Xn Xe) 
wheres 2 
Xe = an Ki (67) 
L= 
In case there exists an odd number of streamtubes, the 
final integral is calculated by: 
E(x) =a +bx o) 
A=Eo | (69a) 
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b= E,-—Eo 


aX, (69b) 
R+ \ 
_ 2 2 
{Eo Ax = A(%24\- X te) + 2 (ee, as Xe ) (70) 


8, The intersection of the characteristic with each stream- 
line, excluding hub and tip, is connected until adjacent 
streamlines divide the flow rate into M-1l equal increments, 
Calculations progress from streamline M-l to streamline 2, 


Bach calculation iterates the increment Ox until: 
K 


Q ()= 27T E(X+ Sx) (dX+d5%) (71) 
le 
Kora | Ecos Ds 
The variables in E(x) are corrected by 5x, 
JK 
Ky = Ri + ($F )i Ox 
= Ri +d5x cos (%— Aid - 
Wrarem Wr + (am 
_ a (Sy }o* 
_ . §X dl Wi (73) 
= Wy +E, 2 ax L 


. =. Xe \GOns 
We= G+ (M4, - G) — 
K indicates the new point of intersection, and i, the 
original point. Eq. (71) is solved by introducing the 
proper limits into Eqs. (66) or (70). 
9, Preliminary calculations for the blade to blade solu- 
tion are included in the computer program for the merid- 


ional plane analysis, Variables of Eq. (34) are grouped, 
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and the groupings computed at the corrected intersections 


of the characteristic and streamlines, 


DW noes > ( 20 at’) cos (> (75a) 


2 
DWeune = BWmtana+ oR (75b) 


C, Computer Programs; ROTOR 1 and LEDGE 1 

1. The preceding methods are translated into FORTRAN com- 
puter language compatible with the CDC 1604 digital computer, 
Program ROTOR 1 calculates a complete solution for one 
characteristic, from the location of the characteristic 


to the correction of assumed streamlines and the relative 


velocity profile, The solution procedes in an orderly 





fashion, much like the development of the preceding section, 
Subroutines are used for repetitive calculations, Control 
of decisions, iterations, and progressive development is 
maintained in the main program, Subroutine PSTAR locates 


the characteristic, accounting for all signs of dN*, Sub- 


ae | 
am é 


by finite differences, Subroutine ANGLE calculates (3 and 


routine DDM calculates the first derivative of ©, 


a Subroutine COY calculates the elements of the coef- 
ficients, Y, except K(n). The coefficients, Y, are calcu- 
lated in the main program, The relative velocities are 
Calculated in subroutine RELVEL, and the flow rate in sub- 
routine FLOW, Control of the iterations is maintained in 
the main program, Iterations for streamline corrections 


are made in the main program, using subroutine FLOW, 
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2. Entrance conditions are calculated in a separate pro- 
gram, since these computations are not required after the 
locations of streamlines on the leading edge are deter- 
mined, Program LEDGE 1 calculates K(n) as defined by Eqs, 
(55), (56), and (57), with the exception that cosX is 
omitted, This variable is a function of position ona 
streamline, and therefore is included in the calculations 
for ¥ 36 The results of LEDGE 1 are used as inputs to 
ROTOR l. 

It should be noted that the first term of K(n) in Eq. 
(12) is omitted, since Ofen O in this analysis, This term 
is introduced as input (DHEDN) into ROTOR 1, It is intend- 
ed that this term would be calculated separately from given 
thermodynamic data at the entrance, 
3, Definitions, flow diagrams, and program listings are 
included in the Appendix, Many control, indexing, and 
grouping names are undefined, Their usage may be inter- 
preted from the developments of the preceding section or 


from the program listings, 


D. Results of ROTOR 1 

This program computed properly, but the results were 
unacceptable, Streamline corrections were large (up to 
1/2 inch), imposing improbable velocity distributions, 
The characteristics appeared to be properly located, how- 
ever, the flow angles, [> , were somewhat arbitrarily 


Gistributed along normals and characteristics, 
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A number of test print-outs were made in order to 
locate the inaccuracies, (The listing of Rotor 1 in the 
Appendix includes test print-out instructions.) It was 
discovered that the derivatives of © , calculated by Eqs, 
(39), (40), (53), and (54) were inaccurate, The second 
derivatives were particularly incorrect, The methods used 
are considered to be sound, however, the degree of accuracy 
primarily depends on the accuracy of the input data, In 
this analysis, the © =constant lines were Originally plot- 
ted from approximations, Errors introduced in plotting 
were compounded by graphical interpolations (and extrapola- 
tions) for 9 values at grid points, The erroneous © de- 
rivatives were introduced in the computations for the co- 


efficients, Y, resulting in an improbable flow solution, 


E. Alternate Meridional Plane Analysis 

1, It was decided to eliminate the derivatives of © from 
the calculations, therefore, it was necessary to introduce 
known values of (= - it is reasonable to assume that these 
data would be available, or could be computed, from detailed 
drawings of the impeller. For this analysis, [> distribu- 
tions along the hub and tip contours were obtained from data 
used in the original design work, However, /3 cannot be 
specified on the internal streamlines until the streamline 
is constructed, A distribution of [> from hub to tip must 
be assumed, 


2, {2 along the hub and tip contours is shown in Fig. 9, 
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It is assumed that (4 varies linearly along normals be- 
tween these contours, 
GO, = Py > (Br -Gu) Os (76) 
Tv 
The function D., of Eq. (49) becomes: 


D> = \ al(R+an p>») 
KR Cos d af 





(77) 


=tan § ba an (> Gand+ B eon/> 
R 2 e085 


The distribution of tan @ between two grid points on a 


streamline or normal is assumed to be linear. 


dian (s_ tan Six - tan So 
Om 


(78a) 
AN 4 
stan (* rr ae 
at the i" grid point: 
dten B _ tan [Sin —ten Pi-r 
OM, 2 DN, a 
At a point, P*: 
atan > ( Stans _ dtan@\ an iS 
oy O Vix oN: AN, 


The angle § is still measured on the © =constant curves, 
Leading edge calculations are also modified, A linear 
distribution of (3 along the leading edge is assumed in the 


same manner as Eq. (76). 


F, Computer Programs: ROTOR 2 and LEDGE 2 


1. Frogram ROTOR 2 is a modification of ROTOR 1, with the 


new calculations for /3 replacing the © derivatives, The 
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order of some solutions has been changed to simplify the 
mechanics of the program, Subroutines DTD and ANGLE are 
eliminated, The values of the angles, ZY , are computed 
in the main program or in subroutine PSTAR, Subroutine 
PSTAR is modified to reflect these changes. Subroutine 
DDL replaces COY, and calculates the derivatives along the 
©=constant curves, Curvatures are calculated in PSTAR, 
Subroutines RELVEL and FLOW are unchanged, Program LEDGE 
2 is a simplification of LEDGE l, reflecting the introduc- 
tion of /3 data, 

Both programs are diagrammed and listed in the Appen- 


dix, Only new or modified variable names are defined, 


G, Results of ROTOR 2 
1. The complete meridional plane analysis was solved by 
ROTOR 2, Nine characteristics were generated, Fig, 4 
is a tracing of the actual construction for the first two 
Characteristics, The grid networks for both solutions are 
indicated, Two modifications are made to the first set of 
estimated streamlines, The first modification reflects the 
changes effected by the calculations for characteristic Cl, 
Each change required recalculations in LEDGE 2, The re- 
sults of the last calculation were held constant throughout 
the remaining solutions, 

Data for each characteristic are listed in Tables II-1l 
through II-9, The "starting normal" originates from the 
"Starting point" for each calculation, The "starting point 


number" indicates the position of the normal in the grid 
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system for each solution, The intersection of a character- 
istic with each streamline is located by measuring the dis- 
tance DNSTAR (dN*) from the intersection of the starting 
normal with the streamline, Streamlines are corrected by 
measuring the distance DELTA (8x) along the characteris- 
tic from the original intersection, P*, Relative velocities 
are listed for the final intersection point, Additional 
data are listed as checks on the calculations, 

For example, the calculations for characteristic C8 
showed a discrepancy between the correction length, 9 x, 
and the computed radius, in locating the new streamlines, 
Test print-outs revealed errors in the calculated grid 
lengths, dM, These errors had accumulated, resulting in 
somewhat radical streamline deviations and questionable 
velocities, dM had been calculated by Eq. (42a), which, 
at the time, was used for angles, rr , of 60 degrees or 
less, This limit was changed to 45 degrees and Eq, (42b) 
applied, The results were satisfactory, 

2. The completed system of computer characteristics and 
streamlines is shown in Fig, 10, True normals are con- 
structed perpendicular to these streamlines, Streamlines 
have been faired and extended to the trailing edge, Data 
points are indicated, 

The distributions of relative velocity along odd num- 
bered streamlines are plotted in Fig, ll, Values at the 
leading edge were hand calculated, Velocities at the trail- 


ing edge were extrapolated on a large scale plot, Curves 


34 








—! ' 


are faired through data points without smoothing, The rela- 
tive velocity profiles on the odd numbered normals of Fig, 
10 are cross plotted in Fig, 12, 

3, Computer run time for one ROTOR 2 solution was approxi- 
mately two minutes, Construction and data preparation time 
was about four hours, Total preparation and computer time 


for LEDGE is approximately two hours, 


H, Blade to Blade Analysis 
1. Most of the preliminary calculations have been accomp- 
lished in the hub to tip solution, The values of W, DW ook! 


and DW are plotted from the results in Tables II, and 


func 
the curves smoothed, Data are read from these curves at 
equal increments along each streamline, Extrapolated data 
are used at the leading and trailing edges, An increment 
of 1/2 inch is used in this solution, 


The derivatives of DW are computed by the five- 


func 
point difference formulas of Eqs. (39). The relative 


velocity difference, AW, is computed by a form of Eq. (34): 
AW = DWeoeg , d(DWeuwe) (81) 
2 am 


The relative velocity distributions on the pressure and 
suction blade surfaces (driving and trailing surfaces ina 
compressor) are calculated by Eqs. (27). 

2, Program BLADE performs these calculations along each 
streamline, from leading edge to trailing edge, Computer 
run time was approximately 30 seconds, Preparation time was 


about six hours, 
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3. Results are compiled in Tables III-1 through III-9, 
The relative velocity distributions along the hub, tip, 
and mean streamlines are plotted in Figs, 13, 14, and 15, 
The linear velocity profiles across the blade channel, at 


various stations along the mean streamline, are shown in 


Fig, 16, 
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IV. DISCUSSION 


A, Theoretical Assumptions 

1, The assumption of axial symmetry is in Keeping with 
sound engineering practice, Although changes in the cir- 
cumferential direction are ignored, incremental pressure 
changes across the blades are considered by introducing 
the blade force, Foe This assumption is also made in 
Refs, 1 and 2, The general theory, postulated by Wu in 
Ref. 4, does not admit to axial symmetry, Wu modifies the 
development of Ref. 2 to allow for deviations from axial 
symmetry, accounting for these deviations with a thickness 
factor, B. This factor is conveniently used as an inte- 
grating parameter in defining a flow function, Bis 
interpreted as being proportional to the thickness of a 
stream sheet containing an arbitrary stream surface, and 
to blade thickness distribution, 

This theory is sound, and would produce accurate re- 
sults if explicitly applied. However, the factor B must 
be assumed, either arbitrarily or from available data. An 
example of an assumed distribution of B is found in Ref, 6, 
A blade to blade analysis is made, using the distribution 
as input, This report states that such data are obtained 
from a meridional plane analysis, presumably by the methods 
of Ref, 4. Yet, Ref, 4 states that B is best calculated 
from the results of a blade to blade analysis, Wu, in Ref, 
13, uses a value of B equal to one, which imposes axial 


symmetry, He states that methods for estimating B were 
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unavailable at the time, This reference also shows that 
the changes effected by this assumption are small. The 
advantage of using approximate values of B, over the 
axisymmetric approach, is doubtful, 

2. The assumption of a steady, inviscid, adiabatic flow 

is generally accepted, Exclusion of any one of these fac- 
tors introduces highly complex theoretical considerations, 
which are difficult, if not impossible to apply. The 
assumption of isentropic flow within the rotor is reason- 
able, since known exceptions are considered, The specific 
exceptions, included in this theory, are the entropy changes 
caused by discontinuities of the leading and trailing edges, 
The analysis of this thesis is based upon the premise that 
these entropy changes are small and can be ignored, This 
is also implied in Ref, 1, which analyzes flows which con- 
form to the prescribed blade angles at the leading and trail- 
ing edges, The analysis of this thesis is not limited by 
this implication, Entrance conditions are not restricted 
to design conditions, and Eq. (25) is generally applicable, 
3, An incompressible solution was developed in this thesis 
in order to simplify the presentation, Compressible solu- 
tions can be solved by this method if the volumetric flow 
rate is replaced by the mass flow rate, Density changes 
along a characteristic can be computed by thermodynamic 
relations which equate the equilibrium conditions at any 
point, P*, to Known or computed entrance conditions, 


4, <A simplified blade to blade solution is used in this 
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thesis, The assumption of a linear velocity distribution 

is not unique, Ref, 1 assumes a linear pressure variation 
across the blade channel, Ref, 6 validates this assump- 
tion for design conditions, The results of Ref, 3 show 
near-linear velocity distributions for both the compressible 


and incompressible solutions, 


B, The Design Analysis 


1, The design analysis was discontinued because of time 
considerations, However, it was shown that the theory 
could be applied to the design solution as well as the 
direct solution, Difficulty arose in prescribing practi- 
cal blade characteristics, One method would be to assume 
various (3 and € distributions as functions of stream- 
line and normal directions; then generate a blade shape 
for each combination, and select the most practical solu- 
tion, Selection would be contingent upon the results of 
a direct analysis, similar to the method developed in this 
thesis, This approach appears to be guite long, 

The question is whether or not this development is 
Warranted, A design method, based on the theoretical de- 
velopment of Ref, 1, is presented in Ref, 5, The hub pro- 
file is calculated from successive flow solutions, This 
approach is modified by Ref. 7 to provide for arbitrary 
entrance and discharge conditions, In effect, the initial 
conditions establish the © =constant curves, The solution 


provided the upper limit of these curves, The design 
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problem is practically solved by the initial conditions, 
All that remains is a system of repetitive calculations, 
which produces an acceptable tip profile, In comparison, 
the method proposed in this thesis would generate blade 
shapes, rather than tip contours, It seems more logical, 
and perhaps more expeditious, to establish feasible physi- 
cal contours and analyze the flow under prescribed condi- 
tions, The contours could then be altered to eliminate 
undesirable flow phenomena, and to produce the desired 
cerformance, This is precisely what is done in the direct 
approach, Therefore, the design, or indirect, method is 


considered to be unnecessary, 


C, Methods of Solution 

1, The developments of Refs, 8 and 9 present a relatively 
simple approach to the direct problem when compared to Ref, 
4, The establishment of characteristic solutions provides 
a method, compatible with high speed machine computations, 
A limited number of applications of Ref, 4 have been pre- 
sented, None have attempted a solution of the complete 
theory, but have reduced the problem to solutions for 
specific applications, Ref. 6 deals with the blade to 
blade analysis, Refs. 5 and 7 are concerned with the d32- 
sign problem, and consider the methods proposed in Ref, 4 
to be impractical for engineering applications, In Ref, 

13 the methods of Ref, 4 are simplified by assuming axial 
symmetry, and an incompressible solution is used as a basis 


for the compressible solution, It is concluded, therefore, 
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that the approach presented in this thesis is theoretically 
sound, 
2. The accuracy of the solution of any theoretical analy- 
sis depends upon the precision of the numerical methods used 
in computations, and in the accuracy of the input data, It 
is difficult to assess the inaccuracies introduced by the 
numerical approximations in this method, It would be im- 
possible to isolate each error, even if test results were 
available for comparison, Gross inaccuracies, such as those 
revealed by ROTOR 1, can be detected, Small errors may can- 
cel, or large errors may ensue from an accumulation of small 
inaccuracies, 

For example, the assumptions of linear distributions 
of curvature and flow angles are questionable in the vicinity 
of the axial midpoint of the test impeller, Hub and tip cur- 
vatures (Fig, 6) follow different patterns in this region, 
The odd distribution of hub curvature is difficult to ex- 
plain, The difference between hub and tip blade angles is 
a maximum in this region, (Fig. 9) These two quantities 
are multiplied in the calculations of the coefficient, Y, 
Difficulties were encountered in obtaining accurate results 
along characteristics in this region, indicating that the 
errors introduced by these assumptions may have multiplied, 
Such reasoning is not conclusive, however, since the results 
of the complete solution do not indicate excessive errors, 

Many interpolations were used in this solution, The 


accuracy of each depends on the grid spacing, The spacing 
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of normals was limited to a maximum of .3 inches, Nine 
streamlines are considered sufficient, A finer network 
would certainly increase the accuracy of each calculation, 
however, the process of measuring data would become tedious, 
The reliability of all esienl atens"aws depends upon 
the accuracy of the measured data, The meridional plane 
was drawn to double scale, Lengths were measured to an ac- 
curacy of about ,05 inches, The angles,  , were measured 
with a drafting machine to an accuracy of 5 minutes, (.1 
degree) This accuracy is reduced to about .3 degrees, since 
measurements were made on curved lines, Interpolations 
were made for & to the nearest tenth of a degree, but an 
accuracy of 1/2 degree is more reasonable, 
3, The general approach to the meridional plane solution 
can be applied to obtain the desired degree of accuracy, 
within the limits outlined in the preceding discussion, A 
finer grid will increase accuracy for a particular charac- 
teristic solution, A solution can be repeated, using the 
previous results to establish new initial conditions. The 
number of characteristic solutions used is arbitrary, 
Therefore, a detailed investigation can be conducted in dis- 
crete regions in a rotor channel, or a complete solution can 
be rapidly calculated, using a few characteristics, 
Theoretically, the method can be used to calculate 
characteristic solutions outside the rotor channel. In this 
case the defining angles of the blade are set equal to zero, 


and the characteristics coincide with normals, In stationary 
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cascades, the rotational effects are eliminated, These 
solutions are excluded in the computer solutions of this 
thesis, because of the schemes adapted in performing nu- 
merical interpolations, Solutions can be calculated in 
close proximity to the leading and trailing edges if fic- 
titious blade characteristics are used outside these 


boundaries, 


D, Results 

1, The results of the first method of solution (ROTOR 1) 
were unacceptable, This was partially due to the numerical 
methods used in calculating the © -derivatives, However, 
inaccuracies of the ©=constant curves are considered to 
be the primary fault, It is believed that accurate © =con- 
stant curves can be constructed from three-view drawings, 
Under this assumption, ROTOR 1 can be utilized with accep- 
table accuracy, 

2. The results of ROTOR 2 are considered valid, Althoush 
test results are not available for comparison, the logic 

of the preceding discussion is sufficient proof, in that: 

a. The theory is sound, 

b, The assumptions are justified for engineering 
applications, 

c, Numerical methods are designed to minimize inac- 
curacies in each calculation, so that errors will 
not accumulate, 

ad, Data are accurately measured, 

3, The computed system of streamlines is regular and follows 
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a logical development through the rotor, Fig, 10 is a 
smooth plot of this system, Deviations from data points 
are small, Maximum deviations occur in the region of rapid 
curvature and blade angle changes, 
4, The acceptable system of streamlines demonstrates the 
consistency of the solution, The accuracy of the method is 
measured by the resultant relative velocity distribution, 
This distribution completely defines the flow, Both pres- 
sure distribution and power delivered, or required, are 
calculated directly from the velocities, Variations in 
relative velocity indicate areas of possible flow separation, 
The accuracy of the relative velocities at the lead- 
ing edge are computed directly from entrance conditions, 
and are accurate, Variations of velocities along selected 
streamlines are shown in Fig, 11, Proper trends are indi- 
cated for a compressor with deleaded type blades, Most of 
the work input is accomplished in the forward section of the 
rotor, diminishing towards the trailing edge, It should be 
noted that the energy level of the fluid is not affected 
by velocity alone. Eq. (6), together with Eq. (9b), shows 
that the local enthalpy is also a function of peripheral 
speed, OR, Therefore, the deloaded blade effectively 
reduces the deceleration of the relative velocity at high 
radii, so that the work input, thus blade loading, will not 
become excessive, The irregular nature of the velocity 
Gistributions in Fig, 11 might be attributed to inaccura- 


cies in the calculations; however, similar presentations 
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in Refs, 5, 6, and 13 support the existence of these irreg- 
ularities, In contrast, the velocity profiles along nor- 
mals in Fig, 12 are guite regular, 

5. The applicability of these results can be illustrated 
by an investigation of flow separation along a contour, 
Relative velocity changes along the hub, near the leading 
edge, appear to be conducive to separation, The rapid de- 
crease in velocity implies an adverse pressure gradient, 

A separation parameter, K, is derived in Ref, 14, A value 


of K greater than .045 indicates probable separation, 


c\ 
k= Aah at (82) 
e W? R,/5 
Reynold's Number is defined as: 
Re= WS 
VY 


Lets: 
mM 


cos (3 


where: 1 = length along blade from stagnation point 


eS 


Along a streamline, from Eqs. (6) and (9b): 


average blade angle 


. WW K 
Ge So NN pen ts oR a (83) 
For incompressible flow: a ag 
ao aN aR ee 
_wi(-W am ~ wR Fiabe (Vcos@ ) (84) 
KK = Ww s-t 


At a point on the hub, m = 2 inches: 


W = 106 £t, /sec, 
=~ Peele £t./sec,—rLe. 
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= 49° from Fig. 9 
Therefore: K = ,0138 
Thus, flow at this particular point, which appears to be 
in one of the more critical regions for separation, should 
not separate, The relatively high value of K does indicate 
possible separation farther along the hub; however, the 
Change of velocity decreases beyond this point, The ef- 
fects of peripheral speed and curvature on separation are 
illustrated in Eq. (84). The value of K is approximately 
~9092 at the point of maximum curvature, at a distance of 
about .4 inches along the tip streamline, The velocity 
gradient at this point is practically Zero, Although the 
value of K is low, it can be surmised that a nominal nega- 
tive velocity gradient might induce separation, 
6, The results of the blade to blade analysis are exten- 
sions of the more exact meridional plane analysis, The 
relative velocity distributions along streamlines (Figs, 
13, 14, and 15) can also be analyzed for local phenomena, 
which may be aggrevated by the correction, AW, This veloc- 
ity difference is used to calculate the fluid forces on the 
blade with the relations of Section II, A presentation of 
the type shown in Fig. 16 is adequate for these calculations, 
An attempt was made to indicate the correlation between 


the trailing edge velocities in Fig, 14, in order to 
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establish some initial conditions for a wake analysis. 

The dashed line indicates a possible velocity distribu- 
tion, but this is only a guess, Figs, 13 and 15 do not 
demonstrate this convergence, In addition, the velocities 
at the trailing edge are calculated from extrapolated data, 


Therefore, this type of correlation is inconclusive. 


&, Extensions of the Method 
1. The meridional plane analysis sufficiently describes 
the flow for an incompressible analysis, The results might 
be extended to include the calculations of: 

a. Pressure distribution 

b, Power input or power required 

Cc, Separation parameter 

qd, Cavitation parameter 

e. Blade forces 
Most of these calculations depend upon the results of sev- 
eral characteristic solutions, and would be included in 
a separate computer program, Many preliminary calcula- 
tions of point functions (pressure, enthalpy, etc.) can be 
included in the main solution, 
2. The methods should be modified to include character- 
istic calculations outside the rotor passage, This would 
permit continuous solutions through multi-stage machines 
and in unbladed passages, Major modifications would be 
required, since a number of decisions would be needed to 
insure the use of proper blade properties in the vicinity 
of blade boundaries, Input requirements would necessarily 
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become more complex, 

3, The theory and methods could be modified to use non- 
dimensional variables, This would be particularly applic- 
able to design studies, where non-dimensional results of 
preliminary rotor configurations could be directly compared, 
4, This method should be extended to include the compres- 
sible solution, A method is outlined in Ref. 9, 

5. The blade to blade method used in this thesis is prac- 
tical in its simplicity, and is sufficiently accurate when 
used in conjunction with the meridional plane analysis, 

It would be interesting to compare the results of this solu- 
tion with a more exact approach, such as the method of Ref, 
6, The extensions to the meridional plane analysis should 


also be included in the blade to blade solution, 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. It is concluded that: 

1. Theoretical assumptions are based on practical 
consideration, and do not impose excessive limitations on 
the solution, 

2. The meridional plane analysis is developed from 
sound theoretical derivations, which provide a relatively 
rapid and accurate method of solution, 

3, The methods which are applied in solving this 
problem are good approximations, The degree of accuracy 
of each characteristic solution depends on the size of the 
coordinate grid and on the accuracy of data measurements, 

4, The progressive method employed in the meridional 
plane solution can be applied at arbitrary intervals to 
obtain the desired accuracy, 

5, The accuracy of the results, based primarily upon 
theoretical considerations, is sufficient for engineering 
purposes, 

6. Results are not sufficiently accurate to compute 
detailed analyses at the trailing edge, 

7. The results of this method are sufficient to com- 
pletely describe the flow within the rotor, and predict the 
performance of the machine, 

8, This solution is restricted to flows within the 
rotor passage, 

9. The computer program, ROTOR 1, can be used if 


©O=constant curves are accurately specified, 
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10, The use of a direct flow analysis is prefer- 
rable to the inverse method in design studies, 

11, The blade to blade analysis is limited by the 
assumption of a linear velocity distribution, but the 


results are acceptable for engineering purposes, 


Be It is recommended that: 

1. The method be non-dimensionalized, 

2. Computations be added to describe pressure dis- 
tribution, blade force, power, and critical flow phenomena. 

3. The solution be modified to include compressible 
flows, and calculations beyond rotor boundaries, 

4, A more accurate blade to blade analysis be used 


in conjunction with this meridional plane solution, 
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COMPUTER PROGRAMS 
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VARIABLE NAMES 


Name Equivalent to 
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DTDM d 9 /dm 
F F(L) 
MO Mo 
OMEGA W 

RL Ro 
THETA © 
ATHETA © 

VM Ven 


Additions for LEDGE 2 
BETA 3 


ABETA (D 
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Units 
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rad, 


in, 
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in), 

deg. 

rad, 
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deg, 


rad, 
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PROGRAMS ROTOR 1 AND ROTOR 2 





VARIABLE NAMES 


Name Equivalent to Units 
BETA @at P* deg. 
PBETA A at P* rad, 
TBETA tan ( 

XBETA (Sat (P*+ Sx) rad, 
ZBETA @at (P*+ & x) deg, 
BLNO N 

CURV Kn in,~? 
PCURV k, at P* in,7? 
Dl Dy 

D2 D, 

D2TDM2 a*@ /an* 

D2TDNM a? © /anam 

DELN | AN in, 
DELTA 5 deg, 
ADELTA & rad, 
DDELTA S at P* deg, 
PDSLTA Sat Px rad, 
YDELTA Sat (P*+ §x) rad, 
DHEDN OH/ on, 

DM dM in, 
DNSTAR an* in, 
DTDM, DT1, DT2 a © /dam 

DWCO DWooer. 


[Ou 


Name 


DWFUNC 


DA 


GAMMA 
GAMZX 
PGAMMA 
TGAMMA 
AGAMMA 
ALAM 
DLAM 
PLAN 
ALAM 
OMEGA 


QS 


AR 

RA 
THETA 
ATHETA 
THICK 


TKEQV 


Equivalent to 


DWeunc 
ax 
6N/2 
a(x) 
Sat P* 
Bat (P¥*+ 5x) 
at Pe 
tan & 
Wat (P*+ Fx) 
~ 
N 
A at P* 
IN at (P*+ 8x) 
OO 
Q (actual) 
Q (calculated) 
AQ 
§ B(x) ax 
R 
R at P* 
R at (P*+ 5x) 
© 
oO 
At 
At! 
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cu, 


cu, 


cu, 


cu, 


Units 


in, 


in, 


deg, 
deg, 


rad, 


per sec, 
£t,/sec, 
ft, /sea, 
ft, /sec, 
ft, /sec, 
in, 

ing 

in, 

deg, 
rad, 

in, 


in, 





WX 
WSQ 


DWSQ 


ADEL 


Yl, Y2, Y¥3 


N, NN, etc, = 
M, MM, etc, = 
MO = 
NO = 


NS = 


Equivalent to 


W 
dw 
* 
ie at P 
wa at (P*+ 6x) 
2 + 
ve at P 


a? /ax 


INDEX NAMES 


Normal number 


Streamline number 


Units 


ft, /sec, 
ft, /sec, 
ft. /sec, 
ft.,/sec, 


ft. °/sec, * 


ft. °/sec, *-in. 
in, 


in, 


in, 


Total number of streamlines 


Total number of normals 


Starting normal 


Changes and additions for ROTOR 2 


ABSTA 
APBETA 
DTBN 


DTBM 


(3 


(2 at P* 
a(tan (> )/dn 
d(tan(? )/am 


lO6 


rad, 


deg. 


— 7 


— A 
ome 


Sa 





S.R. 


S oR. 


S.R. 
for 


Sik. 
for 


Program ROTOR |! 


Main Program 


Read input 


Angles to Radians 


Calculations at Initial Point: 








OTTO < 
ANGLE 

Calculate OM 

Calculate odN* 
PSTAR S.R. PSTAR 
+ ON® for -dN* 
OTO oo eG 
ANGLE le 

dN* = O 

PSTAR <a S.R. PSTAR ° 
+x! far -y! 





> 


o 





a 










Print Results 
at P* 


Yes 


S.R. COY 


& 





Calculate Y's 


SR. RELVEL 


Relative Veloetty 
Distribution 
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Calculations 
at 


. Q 
Pr and Pe 


Inittalize to Correct 
Streamlines 


QAF = QA (M ./M,) : 





Calculate QC 
Re ee 
S.Re FL i) 


Qe———————-—(-) 0)__-»-(4) 
(+) 


8x = se 


S.R. FLOW Calculate QC 


: 






SiR. FLOW A er Calculate QC 


e's 


Sp 





oa 


Calculations for 
Blade to Blade 
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